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ABSTRACT 

We present a new determination of the local volume-averaged star formation rate 
from the 1.4 GHz luminosity function of star forming galaxies. Our sample, taken 
from the B < 12 Revised Shapley-Ames catalogue (231 normal spiral galaxies over 
effective area 7.1 sr) has ~ 100% complete radio detections and is insensitive to dust 
obscuration and cirrus contamination. After removal of known active galaxies, the 
best-fit Schechter function has a faint-end slope of —1.27 ±0.07 in agreement with the 
local Ha luminosity function, characteristic luminosity L* = (2.6±0.7) x 10 22 W Hz -1 
and density </>* = (4.8±1.1) x 10 -4 Mpc -3 . The inferred local radio luminosity density 
of (1.73 ± 0.37 ± 0.03) x 10 19 W Hz -1 Mpc -3 (Poisson noise, large scale structure 
fluctuations) implies a volume averaged star formation rate ~ 2x larger than the 
Gallego et al. Ha estimate, i.e. pi.4GHz = (2.10 ± 0.45 ± 0.04) x 10 -2 M 0 yr -1 Mpc -3 
for a Salpeter initial mass function from 0.1 — 125 M 0 and Hubble constant of 50 km s -1 
Mpc -1 . We demonstrate that the Balmer decrement is a highly unreliable extinction 
estimator, and argue that optical-UV SFRs are easily underestimated, particularly at 
high redshift. 
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1 INTRODUCTION 

Some of the most ambitious and widely-cited extragalactic 
observations in recent years have been the constraints made 
on the star formation history of the Universe. Madau et al. 
(1996) used the ultraviolet luminosity density as a tracer of 
the comoving star formation rate (SFR), and by combining 

2 > 3 Lyman dropout surveys (Steidel et al. 1995) with the 
0.1 < 2 < 1 Canada-France Redshift Survey (CFRS; Lilly 
et al. 1996) they inferred a 2 ~ 1.5 peak and subsequent 
redshift cut-off in the SFR. If correct this is a major result, 
implying an empirical determination of the dominant epoch 
of metal production in the Universe (also suggestively close 
to the peak in QSO number density) which can be com¬ 
pared directly with hierarchical models of galaxy formation 
(Baugh et al. 1997). The dereddened CFRS ultraviolet SFR 
quoted in Lilly et al. (1996) appears to asymptote to the 
local (dereddened) Ha estimate from Gallego et al. (1995). 
However ultraviolet estimates are very sensitive to dust ob¬ 
scuration, with reddening corrections ranging from factors 
of order 2 (Madau et al. 1996) to 10 (Meurer et al. 1999). 
Indeed there have been several claims that less reddening- 
dependent measures yield substantially higher star forma¬ 
tion rates: for example, sub-mm galaxy surveys both behind 


lensing clusters (e.g. Blain et al 1999a,b; Ivison et al. 2000) 
and in the field (e.g. Hughes et al 1998, Barger et al. 1999a,b, 
Lilly et al. 1999a,b, Eales et al. 1999, 2000, Peacock et al. 
2000, Serjeant et al. 2001; Fox et al. 2001 in preparation, 
Scott et al. 2001 in preparation) find star formation rates 
at 2 > 3 comparable to dereddened UV-selected samples, 
but the galaxy populations which appear to comprise the 
sub-mm point sources overlap little with the UV-selected 
star forming population; ISO observations of the Hubble 
Deep Fields at 7/rm and 15/tm (Rowan-Robinson et al. 1996, 
Oliver et al 2001 in preparation) find systematically higher 
star formation rates than inferred at 2 ~ 0.5 — 1 in the 
optical; the Ha luminosity density in the CFRS (Tresse & 
Maddox 1998) implies a z ~ 0.3 star formation rate 2 — 3 
times higher than that inferred from de-reddened ultravio¬ 
let estimates; and the evolving radio luminosity density is 
also consistent with comparable amounts of star formation 
from UV samples and inaccessible to UV samples (Haarsma, 
Partridge, Windhorst & Richards 2000). 

In this paper we will show that the local Ha star forma¬ 
tion rate is itself also a significant underestimate, by compar¬ 
ing it with the local radio luminosity density of star forming 
galaxies. This paper is structured as follows. In section we 
present the sample used in this paper, and corresponding se- 
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Figure 1. Local 1.4 GHz radio luminosity function of Re¬ 
vised Shapley-Ames spiral galaxies, with known AGN from 
Veron-Cetty & Veron (1993) catalogue excluded. Errors are the 
si/v max ± V S1 /V3U range in the case of > 1 galaxy in the 
bin; for bins with a single galaxy the errors represent ±1 a lim¬ 
its on the appropriate Poisson distribution mean. The full line 
shows the Schechter function fit described in the text, and the 
short dashed line shows the original Condon (1989) luminosity 
function. The long dashed line shows the best-fit Ho luminosity 
function of Gallego et al. (1995), scaled using equations j| and 

lection criteria. Section presents a calculation of the local 
radio luminosity function and luminosity density. Implica¬ 
tions of these results are discussed in section [I| We assume 
a Hubble constant of Ho = 50 km s -1 Mpc -1 , qo = 0.5. 


2 SAMPLE DEFINITION 

Unlike the far-IR, Ha and UV, the radio luminosity den¬ 
sity traces both obscured and unobscured star formation, 
as well as being free of cirrus contamination which affects 
far-IR estimates at faint luminosities (Condon, Anderson & 
Helou 1991). At low radio frequencies (~ 1GHz) the lumi¬ 
nosity is dominated by non-thermal emission, plausibly from 
cosmic ray electrons accelerated by supernovae, which can 
be calibrated to the recent star formation rate using the su¬ 
pernova rate in the Galaxy (Condon 1992; section js|). The 
low frequency radio luminosity density is therefore an ideal 
estimator of the local volume-averaged star formation rate. 

The B < 12 spiral galaxies in the Revised Shapley- 
Ames catalogue were mapped at 1.4GHz by Condon (1987) 
with the VLA to limiting flux densities of (typically) 0.1 —0.5 
mjy, obtaining detections of virtually 100% of the sample. 
The survey spans a sufficiently large cosmological volume 
for large scale structure density fluctuations to be negli- 


Table 1. Local 1.4GHz luminosity function parameters: 
<p(L)dL = </>»(L/L*)“ exp(—L/L*)d(I//L*) Mpc -3 . The large 
scale structure density fluctuations are estimated to be around 
2% on the scale of the survey, and are not incorporated into the 
tabulated errors. We assume a Hubble constant of H o = 50 km 
s -1 Mpc -1 . 


Parameter 

min-x 2 value 

drier range 

0* (Mpc -3 ) 

4.9 X 10 -4 

(3.7 - 5.9) x 10 -4 

L* (W Hz -1 ) 

2.8 x 10 22 

(1.9 - 3.3) x 10 22 

—a. 

1.29 

1.21 - 1.35 

L*</>*r(2 + ot) 

(W Hz -1 Mpc -3 ) 

1.73 x 10 19 

(1.35 - 2.08) x 10 19 


gible in comparison with the shot noise on the luminosity 
density and overall number density (Condon 1989; Oliver, 
Gruppioni & Serjeant 1998): using the Peacock & Dodds 
(1994) power spectrum we estimate density fluctuations of 
less than 2% on the scale of the survey, which agrees with 
the uniformity of the < V/Vjnax > statistic derived by Con¬ 
don (1989). Nevertheless, as a precaution the Local Group is 
excluded by eliminating the 10 galaxies within 1.7 Mpc, and 
the Virgo Cluster excluded by removing the 27 within 10° 
of the cluster centre. These exclusions leave a sample of 267 
spiral galaxies. With the correction for galactic extinction, 
the resulting survey has an effective areal coverage of ~ 7.1 
sr (Condon 1989). 

AGN contribution could potentially introduce a sig¬ 
nificant systematic into estimates of the radio luminosity 
density; indeed, Ho et al. (1997) have shown that at least 
50% of spiral galaxies show some evidence for activity when 
a very low detection level is used and the stellar absorp¬ 
tion from the host galaxy is properly corrected. However 
the appropriate question is not whether an active nuclei are 
present, but whether the nuclei contribute significantly to 
the low-frequency radio fluxes. Condon, Anderson & Helou 
(1991) found that removing the Seyfert galaxies fisted in 
the Veron-Cetty & Veron (1993) catalogue from the IRAS 
Bright Galaxy Sample and RSA samples, also eliminated 
the outliers in the radio-far-IR correlation. Although such 
AGN segregation is somewhat arbitrary, this is strong cir¬ 
cumstantial evidence that the remaining galaxies represent 
a more homogeneous population. We can reasonably inter¬ 
pret the remaining galaxies to be those in which the radio 
and far-IR are dominated by star formation. We therefore 
eliminated the 36 galaxies fisted as Seyferts in the Veron- 
Cetty & Veron (1993) catalogue, differing somewhat from 
the approach of Condon (1989) (with a non-negligible ef¬ 
fect around the break of the luminosity function; section 
[j[). The exclusion of AGN neglects any contributions from 
circumnuclear star formation in these objects, and any star 
formation in ellipticals is also explicitly excluded, so the lu¬ 
minosity density from star formation derived in this paper 
should strictly speaking be treated as a lower limit. Fur¬ 
thermore, it is important to note that none of the SFR from 
UV estimates from the CFRS, Lyman drop-out samples and 
the Ha luminosity density of Tresse & Maddox (1998) (and 
indeed most SFR studies) excludes any narrow-line AGN 
contribution. 






The local star formation rate and radio luminosity density 3 


3 THE RADIO LUMINOSITY FUNCTION 

In figure [I] we plot the 1/Unax estimate of the radio lumi¬ 
nosity function (Schmidt 1968) for this sample. Accessible 
volumes were limited on one side by the Local Group red- 
shift limit (section ||), and on the other by the maximum 
redshift at which an otherwise identical galaxy would be ob¬ 
servable with these flux limits. (In practice the optical limit 
dominates the selection.) As noted in Oliver, Gruppioni & 
Serjeant 1998, the presence of two flux limits (optical and 
radio) does not affect the 1/Unax estimate of the luminos¬ 
ity function, as long as (a) the minimum z max redshift of 
the two flux limits is used, and (b) the luminosity function 
is corrected for any parts of the (radio, optical) luminosity 
plane which are inaccessible at all redshifts. In this case the 
latter consideration is negligible. Bins with only one galaxy 
have errors corresponding to la limits on th e mean of the 
appropriate Poisson distribution; errors are yzll/Ulax oth¬ 
erwise. To avoid two consecutive single-membered bins the 
two faintest bins were combined to make a single bin at 
igi8±o.5^-2 We obtained a best-fit solution with a 

reduced yf of 0.8, tabulated in table 1. By comparison, the 
original Condon (1989) luminosity function has y v = 5.5 
with this data set, due to the inclusion of several AGN in 
their fit. Errors on each parameter were obtained treating 
the other parameters as fixed. 

We have already argued that the </* is insensitive to 
density fluctuations in this sample, though it remains to be 
shown that L*, a and the luminosity density are unaffected. 
In order to demonstrate this, we compared the observed ra¬ 
dio luminosity distribution with that expected for the model 
luminosity function, using a fit to the radio-optical correla¬ 
tion given in Condon (1987). This has the advantage of be¬ 
ing independent of the local density distribution. Using the 
Kolmogorov-Smirnov statistic on the luminosity histograms 
to define a likelihood statistic, we find that our min-\ 2 pa¬ 
rameters are an excellent fit to the luminosity distribution. 

The faint-end slope is in good agreement with that ob¬ 
tained for the local Ha luminosity function by Gallego et al. 
(1995). Note that we should not necessarily expect the slope 
of the radio luminosity function, which is a measure of the 
recent star formation rate, to match that of e.g the B-band 
luminosity function, due to the latter having contributions 
from less recent star formation. It is also important to note 
that our constraint on the faint-end slope is dominated by 
objects with luminosities (3 x 10 -3 ~ 0.3) x L*, and is af¬ 
fected little by the exclusion of fainter bins. The data do not 
fit a flat slope throughout (3 x 10 -3 ~ 0.3) x L* (though 
in the faintest bins a flattening cannot be excluded). This 
“useful” faint-end luminosity range is ~ ldex deeper rela¬ 
tive to L* than that of Gallego et al. 1995, and comparable 
to e.g. that of Tresse & Maddox 1998. The total luminos¬ 
ity density of a Schechter luminosity function is given by 
</*L*r(2 + a), but the errors in these parameters are highly 
correlated. We therefore explored the parameter space to 
find the regions where Pr(\ 2 ) < 0.68, and hence determined 
the ±ler luminosity density range, also listed in table 1. 

At these frequencies the radio luminosity is dominated 
by non-thermal emission from relativistic electrons (Condon 
1992) accelerated by supernovae from massive (M > 5 Mq) 
stars. The lifetimes of such stars are only ~ 10 7 yr, so the low 
frequency radio emission should be proportional to the re¬ 


cent star formation rate. Following Condon (1992) and Con¬ 
don & Yin (1990), this can be calibrated using the Galactic 
non-thermal radio luminosity, obtaining 

SFR( 100 > M/Mq > 5) = M Q yr~\ 1) 

where v is the frequency and a is the non-thermal radio 
spectral index (a = — dlogS^/dlog v ~ 0.8). Assuming a 
Salpeter initial mass function (IMF; ip(M ) oc M -2 ' 35 ) we 
can correct to 0.1 — 125 Mq if we divide by (100 -0 ' 35 — 
5“ 0 ' 35 )/(125 - 0 ' 35 — 0.1 -0 ' 35 ) = 0.18, and can correct from 
1 GHz to 1.4 GHz with the factor 1.4 -0 ' 8 . The non-thermal 
contribution at 1 GHz is ~ 90% of the total luminosity (Con¬ 
don 1992), so we therefore obtain for the star formation rate 

SFR(125 > M/M 0 > 0.1) = Ll ' 4 8.2x [ 10^ Z ^ 

where L 1.4 ghz is the total monochromatic luminosity at 
that wavelength. Our observed luminosity density (table 1) 
therefore corresponds to a volume-averaged SFR of (2.10 ± 
0.45) x 10 -2 Mq yr -1 Mpc -3 . The SFR is dominated by 
objects around L*, peaking at ~ 0.7L* corresponding to 
around 25M 0 yr -1 . 

We can compare this estimate with that obtained from 
the local Ha luminosity density quoted by Gallego et al. 
(1995). For our adopted IMF, Madau et al. (1998) quote a 
conversion of 

SFR( 125 > M/Mq > 0.1) = M 0 j/r _1 (3) 

for unreddened (i.e. dust-free) Ha luminosity densities. Gal¬ 
lego et al. dereddened their Ha fluxes using Balmer decre¬ 
ments, and their resulting luminosity density corresponds 
to a volume-averaged SFR of (8.8 ± 0.1) x 10 -3 Mq yr -1 
Mpc -3 , about a factor x2 lower than our radio estimate. 
We will discuss the reasons for the Balmer decrement failing 
in section a similar systematic offset was seen in the (al¬ 
beit inhomogeneous) samples of Cram et al. (1998), whose 
Ha-derived star formation rates were a few to 10 times lower 
than radio estimates, with the more luminous systems show¬ 
ing the larger offset. Indeed our luminosity function is in 
good agreement with that of Gallego et al. at the faint end 
(figure jl), but shows a consistently higher number density 
around and above our L* where the contribution to the lu¬ 
minosity density is significant. 


4 DISCUSSION 

We have found that even the de-reddened Ha luminosity 
density underestimates the local volume-averaged star for¬ 
mation rate by a factor ~ x2. (Strictly speaking our radio 
SFR is a lower limit, as it excludes the SFR from irregu¬ 
lars and ellipticals.) However, this does not imply that the 
global SFR is dominated by extreme, highly obscured ultra- 
luminous galaxies - indeed we have found the local lumi¬ 
nosity density dominated by only ~ 25 Mq yr -1 systems. 
Neither does it imply, for instance, that spiral disks are pre¬ 
dominantly opaque. It does however require that a signifi¬ 
cant fraction of star formation occurs in the cores of giant 
molecular clouds, which must survive destruction from pho¬ 
todissociation and supernova-driven shocks long enough to 
obscure a significant fraction of the O and B stars. 
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At faint radio luminosities, where we have suggested the 
lowest extinction, the radio-far-IR correlation becomes non¬ 
linear. This is usually attributed to the increasing cirrus con¬ 
tribution in the far-IR. This contribution can be estimated 
using the B magnitudes (Devereux & Eales 1989) but at the 
cost of a greatly increased scatter in the corrected correla¬ 
tion. However, Condon, Anderson & Helou (1991) showed 
that this can be avoided by allowing the cirrus contribu¬ 
tion to also depend on the radio luminosity - in essence, 
the cirrus (of whatever temperature) is partly heated by un¬ 
obscured OB stars. This is thus in perfect accord with our 
lower extinction at the faint end of the luminosity function. 

Conversely, the extinction corrections are significant at 
around L*. One naive interpretation of our high radio SFRs 
would be to invoke a second, high extinction component to 
the SFR. However, two components are not necessary. Con¬ 
sider a simple constant-density model where the dust in the 
star forming regions is well mixed with the Ra-emitting gas. 
Following Kroker et al. 1996 and Thronson et ah 1990, the 
observed Ha flux S of the dusty star forming region is re¬ 
lated to the “intrinsic” Ha flux So which would be observed 
in the absence of dust, by 

S = S 0 {l-e~ T )/ T (4) 


where r ~ 0.7 Ay is the Ha optical depth to the rear of 
the star forming region. Assuming an extinction correction 
of > 1 magnitude was applied to the Gallego et al. (1995) 
SFR (Kennicut 1983, Tresse & Maddox 1998), our radio SFR 
implies 


S > 0.88 x lCT 04 
So ~ 2.1 


0.167 


( 5 ) 


and hence Av it 9. (It is important to stress that this is 
the extinction to the rear of the cloud, and not a “typical” 
extinction.) The observed H/3 flux can also be derived us¬ 
ing t(H/J) ~ 1.45r(Ha). Remarkably, the predicted Balmer 
decrement in this model would imply only Av — 1.1 in this 
model if one (wrongly) assumed a simple dust screen, regard¬ 
less of the true Ay, due to the low-Ay regions dominating 
the observed Balmer line fluxes. This is in excellent agree¬ 
ment with the typical observed Balmer decrements in Ha 
surveys (Kennicut 1983, Tresse & Maddox 1998). The red¬ 
dening corrections derived from optical-UV colours will also 
be underestimates, since the observed spectra are also dom¬ 
inated by the regions with the smallest obscuration. This 
effect will be strongest at high redsliifts, where the rest- 
frame wavelengths are the shortest. Naive reddening cor¬ 
rections could therefore easily introduce redshift-dependent 
biases into the comoving SFR estimates. Even “effective” 
reddening corrections based on typical starburst optical-UV 
SEDs (e.g. Heckman et al. 1998) could miss the most highly 
obscured components unless such SFRs are explicitly nor¬ 
malised to a more isotropic indicator (e.g. Meurer, Heckman 
& Calzetti 1998, Cram et al. 1998). 

Nevertheless, the prospects for reddening-independent 
constraints on the comoving SFR at high -2 are excellent, 
due mainly to the fact that the strong far-IR-radio cor¬ 
relation implies the far-IR can also provide unbiased SFR 
estimates. For example, optical spectroscopy of new sur¬ 
veys conducted with the Infrared Space Observatory, such as 
ELAIS (Oliver et al. 2000), ISO-HDF (Serjeant et al. 1997, 
Rowan-Robinson et al. 1997) and others (e.g. Taniguchi et 


al. 1997, Lemonon et al. 1998), will shortly probe the SFR 
to redshifts of 2 ~ 1 — 2. Unlike the dereddened UV and 
Ha estimates, these mid-IR and far-IR selected samples are 
sensitive to star formation in even the most obscured gi¬ 
ant molecular clouds, with in some cases independent con¬ 
sistency checks available from low-frequency sub-mjy radio 
follow-up observations (e.g. Ciliegi et al. 1999, Gruppioni 
et al. 1999). At higher redshifts still, sub-mm sky surveys 
are already providing strong constraints (e.g. Hughes et al. 
1998; Barger et al. 2000). 

Finally, we note that Madau et al. (1998) could repro¬ 
duce the K-band galaxy counts in an Q = 1, A = 0 cosmol¬ 
ogy. Our results are therefore suggestive of a non-zero A, as 
already suggested by many other recent approaches (see e.g. 
Peacock et al. 2001 and refs, therein). 
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